How does the relaxation of a supercooled liquid depend on its microscopic dynamics? 
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Using molecular dynamics computer simulations we investigate how the relaxation dynamics of a 
simple supercooled liquid with Newtonian dynamics differs from the one with a stochastic dynamics. 
We find that, apart from the early /3-relaxation regime, the two dynamics give rise to the same 
relaxation behavior. The increase of the relaxation times of the system upon cooling, the details of 
the a-relaxation, as well as the wave vector dependence of the Edwards-Anderson-parameters are 
independent of the microscopic dynamics. 
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If the logarithm of a transport quantity, such as the vis- 
cosity 77, of a good glass former is plotted versus Tg/T, 
where T is temperature and Tg is the glass transition 
temperature, it becomes obvious that this temperature 
dependence is not universal since some materials show 
essentially an Arrhenius behavior whereas others show 
a strong non- Arrhenius behavior Q. Also more micro- 
scopic dynamical properties, such as the Raman spec- 
trum, depend strongly on the material, in that, i.e., the 
so-called boson peak is much more pronounced in strong 
glass formers than in fragile glass formers ||^. Thus we 
can say that it is well established that the macroscopic 
as well as the microscopic dynamics of supercooled liq- 
uids is not universal at all and must be considered as a 
material specific property. This insight is of course not 
surprising, since the materials differ in their structure, 
the masses of the individual atoms etc. and thus it can 
be expected that these microscopic quantities will give 
rise to a different relaxational and vibrational dynamics. 

What is much less obvious, however, is how the micro- 
scopic dynamics affects the vibrational and relaxational 
dynamics of the system, i.e. whether the relaxational dy- 
namics is different if the microscopic dynamics is, e.g., a 
Newtonian one or a Brownian one. The answer to this 
question is most important since it will allow to gain 
insight to understand which aspects of the relaxation be- 
havior are, for a given system, universal and which ones 
are not. This information is in turn relevant for test- 
ing the applicability of theories that attempt to describe 
the slowing down of the system upon cooling, i.e. the 
mechanism for the glass transition. 

In real experiments it is of course difficult to investi- 
gate how the microscopic dynamics affects the dynamics 
of the system at long times, since usually it is not pos- 
sible to change the former without also influencing other 
microscopic quantities like, e.g., the masses of the parti- 
cles or the interaction between the atoms. Experimen- 
tally a Brownian type dynamics can be realized, e.g., by 
colloidal fluids |p[, while atomic liquids have a Newto- 



nian dynamics. However, the structure and interparticle 
forces in these two types of systems are quite distinct 
from each other and thus it is not surprising that the 
two corresponding dynamics are different. For computer 
simulations it is, however, most simple to change, for a 
given system, the dynamics and they are therefore ideally 
suited to investigate such questions. The only investiga- 
tion in that direction we know of is a pioneering study by 
Lowen et al. [Q. In that work the authors compared the 
relaxation dynamics of a polydisperse system of charged 
particles which move according to a Newtonian dynam- 
ics to the one moving according to a Brownian dynamics. 
The outcome of that study was that the relaxation dy- 
namics depends, on the time scale of the /3-relaxation, on 
the microscopic dynamics. However, due to the limited 
length of the runs and the lack of statistics no stringent 
test could be made whether or not the a-relaxation de- 
pends on the microscopic dynamics. If such a dependence 
would also exist in the a-relaxation regime it would be in 
contrast to the prediction of the so-called mode-coupling 
theory (MCT) |,§, according to which the relaxation 
dynamics at long times should be independent of the mi- 
croscopic dynamics. Since there is very good evidence 
that this theory gives a reliable description of the dy- 
namics of supercooled simple liquids and even network 
forming liquids, see and references therein, such a 

disagreement between theory and computer simulation 
would be quite disturbing since it shows that certain fun- 
damental aspects of the theory are not correct. Due to 
the availability of better computers and algorithms it is 
today possible to do simulations which are more than 
hundred times longer than the ones done by Lowen et 
al. Therefore it is possible to investigate the question 
how the long time dynamics of a system depends on the 
microscopic dynamics on a qualitatively completely dif- 
ferent level and in the present paper we report the results 
of such an investigation. 

The system considered is a 80:20 mixture of 1000 
Lennard-Jones particles consisting of two species of par- 
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tides, A and B. All the particles have the same mass m 
and the interaction between two particles of type a, /3 G 
{A,B} is given by V^pir) = AeapHcra^/ry^ - [cTap/rf] 
with EAA = 1-0, aAA = 1-0, eAB = 1-5, ctab = 0.8, 
ess = 0.5, and asB — 0.88, with a cut-ofF radius of 
2.5i7q/3. In the following we will always use reduced time 
units with aAA and caa the unit of length and energy, 
respectively (setting the Boltzmann constant ks equal to 
1.0). Time is measured in units of y/a^j^m/ASeAA- The 
volume of the box is kept constant with a box length 
of 9.4. Two types of microscopic dynamics are investi- 
gated: A Newtonian one and a stochastic one (described 
below). In the Newtonian dynamics (ND) we integrate 
Newton's equation of motion with the velocity form of 
the Verlet algorithm with a time step of 0.02. After equi- 
librating the system in the canonical ensemble we turn 
off the heat bath and start the production run in the 
microcanonical ensemble. For this kind of dynamics the 
temperature dependence of the relaxation behavior has 
already been studied in great detail in the temperature 
range 5.0 > T > 0.466 Due to an improvement of 
the hardware and computer codes it is, however, possible 
today to perform simulations which extend over a time 
range that exceed the ones of the previous investigations 
by a factor of eight, giving now a total of 4 x 10^ time 
steps, and thus allowing us to equilibrate the system at 
lower temperatures than has been possible before and 
hence to perform more accurate tests. Therefore we have 
determined the dynamics of the system also at T = 0.452 
and 0.446. 

The stochastic dynamics (SD) we considered is defined 
as follows: Apart from the deterministic forces that orig- 
inate from the interaction potential given above, each 
particle j is also subject to a gaussian distributed white 
noise force f]j{t) with zero mean, i.e., {ffj{t)) = 0, and 
a damping force which is proportional to the velocity of 
the particle. Thus the equation of motion for particle j 
reads: 

m'^i + - rj\) = + rfj{t) , (1) 

I 

where the damping constant is given from the fluctua- 
tion dissipation theorem, i.e. (rfjit) -ffiit')) = 6kBT(^S{t — 
t'). The value of C was set to 10, which is sufficiently large 
that the presented results for the SD do not depend on ^ 
anymore (apart from a trivial change of the time scale) . 
Equations (|l|) where solved with a Heun algorithm [p^ . 
The time step used was 0.008, which is small enough to 
ensure that the equilibrium properties of the system are 
the same for the ND and SD (this was checked at all 
temperatures). Also for this dynamics the runs at the 
lowest temperature extended over 4 x 10^ time steps. At 
the lowest temperatures the SD is so slow that, within 
the time span of the runs, the system does not equili- 
brate. Therefore we used at each temperature the ND 
to equilibrate the system and used the SD only for the 



(equilibrium) production runs. In order to improve the 
statistics of the results we averaged at each temperature 
over eight independent runs (ND as well as SD). The 
temperatures investigated were 5.0, 4.0, 3.0, 2.0, 1.0, 0.8, 
0.6, 0.55, 0.500, 0.475, 0.466, 0.452, and 0.446. 

One of the simplest quantities to characterize the dy- 
namics of the system is the diffusion constant Da of the 
particles, which we calculated from the mean squared dis- 
placement of a tagged particle. For the ND the temper- 
ature dependence of has been determined before ||] 
and it was shown that, at low temperatures, it is given 
by a power-law, i.e. Da oc {T — Tc)'' , a functional form 
predicted by MCT. (Here Tc > is the so-called critical 
temperature of MCT.) It is, however, desirable to com- 
pare the temperature dependence of Da for the ND with 
the one for the SD without making reference to any the- 
ory and thus we calculated the ratio Da nd/Da,sd and 
show its temperature dependence in Fig. ]^. Note that we 
plot this ratio versus T — Tc, were Tc = 0.435 is the crit- 
ical temperature from the MCT analysis instead of 
T, but we emphasize that this representation of the data 
has nothing to do with MCT but is only a conveniant 
way to expand the abscissa at low temperatures. From 
this figure we see that this ratio shows a noticable tem- 
perature dependence at high temperatures but becomes 
essentially constant (within the error bars) for tempera- 
tures T < 0.8. We note that in the temperature interval 
0.446 < T < 0.8 the diffusion constants change by al- 
most three decades Since the ratio Da,nd/ Da,sd 
stays constant to within about 30% we conclude that the 
temperature dependence of the diffusion constants is in- 
dependent of the microscopic dynamics to within a few 
parts in 10^. A similar result is found for the B particles. 

A further quantity that is very useful to characterize 
the dynamics of the system is the (incoherent) interme- 
diate scattering function Fs{q,t) for wave vector q [p^ . 
For the SD we show the time dependence of this quantity 
in Fig. H for all temperatures investigated (solid lines). 
The value of q is 7.20, the maximum of the static struc- 
ture factor for the A — A correlation ||^. It can be seen 
that on lowering the temperature the relaxation behav- 
ior changes qualitatively in that at high temperatures 
the decay of (g, t) is essentially an exponential (apart 
from the time dependence at very short times) whereas 
it shows a plateau at low temperatures. The time for 
which the correlation function decays to zero increases 
quickly with decreasing temperature (note the logarith- 
mic time axis) indicating the dramatic slowing down of 
the relaxation dynamics of the system upon cooling. This 
strong dependence of the dynamics on temperature has 
also been seen in the case of the ND ||] . We emphasize 
that these curves are all equilibrium curves. 

Also shown in the figure is Fg {q, t) for the ND at three 
different temperatures (dashed curves). We first compare 
the SD and ND in the a-relaxation regime, i.e. the time 
regime where the correlation function decays, at low tem- 
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peratures, from the mentioned plateau. From the figure 
we see that at the highest temperature (T = 5.0) the ND 
gives rise to a relaxation that is about a factor of seven 
faster than the one for the SD. This factor increases upon 
cooling the system and reaches about 20-30 at low tem- 
peratures. Thus from the point of view of the absolute 
values of the a-relaxation times t{T), the two types of 
dynamics are very different. {t{T) can be defined, e.g., 
as the time the correlator takes to decay to e^^ of its 
initial value.) However, if we look at the temperature 
dependence of t(T), we come to a different conclusion. 
In Fig. 1^ we also show the ratio ta.nd /ta,sd and find 
that this ratio becomes independent of T for tempera- 
tures smaller than 1.0. Thus the ratio shows a similar 
dependence on temperature as the one for the diffusion 
constants. 

From Fig. ^ we also see that the shape of the corre- 
lation functions in the a-relaxation regime is indepen- 
dent of the microscopic dynamics. This can be recog- 
nized by the fact that the curve for at T = 0.466 for 
the ND lies almost on top of the curve for the SD for 
T ~ 0.55 and the fact that for both types of dynam- 
ics the time-temperature superposition principle (TTSP) 
holds (see Refs. |^ and |ll|), i.e. that for a given dynam- 
ics the shape of the correlation function does, in the a- 
relaxation regime, not depend on temperature. Therefore 
we conclude that not only the temperature dependence 
of the time scale of the a-relaxation is independent of 
the microscopic dynamics but also those details of the 
(a-)relaxation process that give rise to the stretching of 
the correlation function. The same result is found for 
the coherent intermediate scattering function and other 
values of q. 

We now turn our attention to the /3-relaxation, i.e. the 
relaxation regime in which the correlation function is in 
the vicinity of the plateau. From Fig. ^we recognize that, 
for low temperatures, the way the curves approach the 
plateau depends strongly on the microscopic dynamics 
in that for the SD this approach is very gentle whereas 
for the ND it is quite abrupt. This difference does not 
exist for the late /3-relaxation, i.e. when the correlation 
functions start to fall under the plateau. This can be 
seen by plotting the relaxation functions versus i/r(T), 
which is done in Fig. ^ for the SD and the ND for the 
lowest temperature investigated. We see that this scaling 
of time makes the two correlation functions fall on top of 
each other in the late /3-relaxation regime. 

MCT predicts that in the (3-relaxation regime and 
asymptotically close to the critical temperature Tc the 
shape of the correlation functions is given by the so-called 
/3-correlator, a function that can be computed within the 
framework of the theory . The form of this /3-correlator 
depends on one parameter, the so-called exponent param- 
eter A, which has been computed for the present system 
to be A = 0.71 Using this theoretical value of A we 

fitted the correlation curves in the /3-regime with the cor- 



responding /3-correlator and the resulting fit is show in 
Fig. H. It has been shown that, for temperatures slightly 
above T^, the time window in which the theoretical curve 
describes the data well can be extended considerably, if 
also the corrections to the /3-correlator are taken into ac- 
count l^],^ . For long times the leading term of these cor- 
rections is given by a power-law with an exponent which 
can be computed from A. The result of a fit to our SD 
data with the /3-correlator and the first correction at long 
times is included in Fig. ^ as well. We recognize that the 
fit with the /3-correlator alone describes the SD data well 
over about 3.5 decades in time. This time window is 
expanded at large times by about half a decade by tak- 
ing into account the corrections to the /3-correlator (time 
window between the two arrows). From the figure we 
see that the theoretical curves fit the ND data only well 
in the late /3-relaxation regime, whereas they do a poor 
job in the early /3-relaxation regime since the asymptotic 
law are completely obscured by the phonons . For the 
SD, however, also the early /3-relaxation is described very 
well by MCT in that the approach to the (quasi) plateau 
is described well by the /3-correlator. The reason for this 
better agreement is likely that in the SD the phonons are 
strongly damped and thus interfere much less with the 
asymptotic laws of the theory on the time scale in which 
the correlators approach this plateau. 

The last quantity we discuss is the q dependence of /c, 
the so-called Edwards- Anderson or nonergodicity param- 
eter. This quantity is related to the height of the plateau 
in the correlation functions and can be measured, e.g., in 
neutron scattering experiments (see, e.g., Ref. (l^). fc 
is one of the fit parameters for the /3-correlator and was 
thus obtained by performing such fits to the incoherent 
and coherent intermediate scattering function at various 
wave-vectors q. In Fig. ^ we show the g-dependence of 
the corresponding nonergodicity parameters, fc^\q) and 
fc{q), for the SD and the ND for the coherent {A - A) 
and incoherent intermediate (A) scattering function. We 
see that for all q considered the two fc for the SD are 
very close to the ones for the ND, thus showing that also 
this quantity does not depend on the microscopic dynam- 
ics. Also included in the figure are the two curves that 
correspond to the g-dependence of the fc as predicted 
by MCT [^. We see that these two theoretical curves 
agree very well with the ones measured in the simula- 
tion. The reader should appreciate that in the calculation 
of these theoretical curves no free fit parameter was in- 
volved, since the only input to the MCT calculations was 
the temperature dependence of the partial structure fac- 
tors, which were determined from the simulation. Thus 
we conclude that the nonergodicity parameters do not 
depend on the microscopic dynamics and can be calcu- 
lated with high precision from MCT. (We also mention 
that a comparison between the g-dependence of the non- 
ergodicity parameter of the ND and the ones predicted 
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from MCT has been made in Ref. using the fdq) [14] B. Frick, B. Farago, and D. Richter, Phys. Rev. Lett. 64, 

and fc^\q) as determined in Ref. However, as we ^^^^ (1990). 

found out in the course of the present work, these q- 

dependencies were affected by a systematic error which 

originated from the neglect of the corrections of the /3- 

correlator Therefore the curves for the ND shown 

in Fig. 3 are not quite the same as the ones shown in 

Ref. ) 

In summary we conclude that at low temperatures the 
a-relaxation dynamics of this system is independent of 
the microscopic dynamics whereas the early /3-relaxation 
dynamics does depend on it. We find that in the f3- 
relaxation regime the relaxation behavior for the SD is 
in very good agreement with the one predicted by MCT 
and that this theory is also able to give a very accurate 
prediction of the q-dependence of the nonergodicity pa- 
rameters. 
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FIG. 1. Temperature dependence of the ratio of the diffu- 
sion constant for the A particles for the ND and SD (circles). 
Squares: Ratio for the a-relaxation time of the incoherent 
intermediate scattering function. 




FIG. 2. Time dependence of the incoherent intermediate 

scattering function for the SD for all temperatures inves- 
tigated (solid lines) and the ND for selected temperatures 
(dashed lines). 
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FIG. 3. Fs{q,t) at the lowest temperature investigated ver- 
sus roscaled time t/T{T) for the SD (bold sohd Una) and the 
ND (bold dashed line). Dashed line: /3-correlator for A = 0.71. 
Thin solid line: /3-correlator + leading corrections at long 
times. 




FIG. 4. Wave-vector dependence of the nonergodicity par 
rameter for the SD and ND (filled and closed symbols, respec- 
tively). The squares and circles correspond to the incoherent 

and coherent intermediate scattering function, respectively. 
The two solid lines are the prediction of MCT for this sys- 
tem. 
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